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Abstract 

It is pointed out that the moments of phase-space particle density 
at freeze-out can be determined from the coincidence probabihties of 
the events observed in multiparticle production. A method to measure 
the coincidence probabihties is described and its vahdity examined. 

PACS numbers 25.75. Gz, 13.65.+i 

Bose-Einstein correlations, interaction region determination. 



Recently we have proposed P a method for measuring the average of 
the particle density in multiparticle phase space. Generalizing the approach 
used by Bertsch for single particle phase space densities, we defined the 
average density for a final state consisting of M particles as 

{Dm) = M jdXdKW^{K,X). (1) 
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Here Vr(K, X) is the M-particle distribution function and the integrations 
are over the 3M components of the position vectors and over the 3M com- 
ponents of the momentum vectors of the particles. Following Bertsch we 
replaced the emission function depending on the four-vectors K and X by a 
function depending on the three- vectors K and X at some representative time 
to which does not need to be specified. This density distribution represents 
the situation at freeze-out. 

In the present paper the method of measurement is generalized to arbi- 
trary integer moments (/) of the phase space density 



As our argument follows closely that of pi , we shall be very brief. For details 
the reader can consult jT]. 

The method uses the coincidence probabilities defined as follows. An M- 
particle final state could be characterized by specifying the momenta of the 
M particles p^, . . . ,Pjvf Then, however, no two observed events would be 
identical. It has been proposed instead 0-13, to introduce momentum bins 
of finite width and to consider two events as identical, if they have the same 
distribution of particles among the bins. Let us denote the total number of 
events by and the number of states which occurred / times by Ni. If a 

state occurs k > I times it contributes to Ni its weight ^ f ^ • Then the 

number of /-fold coincidences is defined by 



C^^(/) is the ratio of the number of sets of / events with equivalent final states 
to the total number of sets of / events. Thus, for large A^, it is the probability 
of an /-fold coincidence. The special case / = 2 had been considered, in a 
different context, by Ma |H]. The coincidence probability depends, of course, 
on the binning. For a given binning it is in principle not difficult to find from 
a data sample the coincidence probabilities. In practice the limiting factor 
is statistics. 

We will now show that for suitably defined bins the effective coincidence 
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probabilities 



Cm{1 + 1) = ( ) (D^m) (4) 



can be approximated by C'}^^{1 + 1)"'^. 

To see that, we first express the measured 1-fold coincidences (JHl) by the 
3M-dimensional distribution of momenta 

w{K) = jdX W(K,X). (5) 

Let us denote the 3M-dimensional momentum bins by j = 1, . . . , J and their 
volumes by uj. Then the probability that an event corresponds to a point in 
bin j is 



Pj = dK w{K) (6) 
and the average density in bin j 

= ^. (7) 

UJj 

The probability that / events chosen at random correspond each to a point 
in bin j is Pj. Therefore, the probability that / events chosen at random are 
the same is 

Next, we observe that, as seen from (Pj). 

(D') =M^J2 [ d^ P^'+^(K,X). (9) 

In order to proceed further, following P, we restrict ourselves to the 
phase-space distributions of the general form 

l^(K, X) = / G{y./L)w{K). (10) 

^The power of 2tt in appears because for M free particles there is one quantum 
state per volume (27r)'^*^ in phase space. 
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Here X/L is a 3M-dimensional vector with components (Xi— Xi)/Lj., . . . , {Zi- 
Zm)/ Lz- The parameters L^, Ly, L^, Xi, . . . , Zm are in general functions of 
K. They could also be different for different kinds of particles. Function 
G{u) satisfies the conditions 

jduG{u) = l] jduuG{u)=0] jduu^G{u) = l, (11) 
where all the integrations are 3M-fold. These relations imply 



dXG{X/L) = {L^LyL,)^- (12) 
(Xi) = Xi,...,(Zm) =^m; (13) 



a\X,) = Ll- a\Y,) = L'; a\Z,) = for z = 1, . . . , M. (14) 



The first condition ensures that iy(K) is the observed 3M-dimensional mo- 
mentum distribution. The other two yield the physical interpretation of the 
parameters Xi,Yi, Zi, L^, Ly, L^. 

Using the Ansatz (jlOj) and the definition one finds 

Cm{1 + 1) = {2ngi^,f''^J2(f[(^-Wf'^] [ dK w^+\K). (15) 

j \m=l / -^"^J 

with {gi+if^' = JduG'+\u). 

Comparing this formula with formula (jH)) it is seen that if the bin sizes 

are 



= 1 iLxLyLz)j 

we obtain 



V, ^ r dK[w(K)y 

Eq. p6p is rather general and can be applied to an entirely arbitrary 
discretization procedure. In the simple (but probably most practical) case 
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when LxLyLz does not depend on K, one can choose bins of constant lengths 
Aj;, Ay, A2 along the axes K^, Ky, K^, respectively. Then the condition for 
the size of the bin is 



A.A,A, = (18) 



{2m 

■'xLyJ-lz 

and Uj = [A^Aj^A^]^. 

These results, being very similar to that of P], invite several comments 
which were elaborated there at length. Here we only repeat the main points. 

As a first approximation we may omit the ratio of averages in the square 
bracket, giving Cm{1) = Cm^(0- This approximation becomes exact when 
function w(K) is constant within each bin j. The condition is on the volume 
of the bin, but does not constrain its shape. This freedom can be used 
to improve the approximation. Bins should be chosen narrow in directions 
where function w(K) changes rapidly and broad in directions with little or no 
variation of w(K). For given variability of w(K) our approximation improves 
when the bins become smaller, i.e. when the product L^LyLz increases. Thus 
we expect the method to work much better for central heavy ion collisions 
than e.g. for e"'"e~annihilations. With very good statistics, or a reliable 
Monte Carlo event generator, one could estimate the term, which here we 
replaced by unity, and thus improve the approximation. In order to find the 
bin size it is necessary to know the volume in coordinate space {L^LyL^) 
and the integral jdu G\u). Information about the volume may be provided 
by interferometric measurements. The integral is not very sensitive to the 
shape of the distribution. For a rectangular box 271 gi = 11/ -\/3. For Gaussians 
27r^i = y2^/[y/]V('-i). 

In conclusion, we have generalized the results of [I] to measurements of 
the higher moments of the phase-space density of particles produced in high- 
energy collisions. It turns out that these moments can be approximated 
by the measured coincidence probabilities of the appropriately discretized 
events. The optimal discretization is suggested and shown to depend criti- 
cally on the volume of the system in configuration space. Thus the actual 
measurements would require this additional information to be effective. 
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